Astron. Nachr. / AN 999, No. 88, 789 -|794](2006) / DOI please set DOI! 



MUCHFUSS - Massive Unseen Companions to Hot Faint Underlumi- 
nous Stars from SDSS 

S. Geier 1 *, V. Schaffenroth 1 , H. Hirsch 1 , A. TiUich 1 , U. Heber 1 , P. F. L. Maxted 2 , R. H. 0stensen 3 , B. 
N. Barlow 4 5 , S. J. O'Toole 6 , T. Kupfer 1 7 , T. Marsh 8 , B. Gansicke 8 , R. Napiwotzki 9 , O. Cordes 10 , S. 
Miiller 1 , L. Classen 1 , E. Ziegerer 1 , and H. Drechsel 1 

j 1 Dr. Karl Remeis-Observatory & ECAP, Astronomical Institute, Friedrich-Alexander University Erlangen-Nuremberg, 

I ■ Sternwartstr. 7, D 96049 Bamberg, Germany 

' 2 Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK 

£SJ | 3 Institute of Astronomy, K.U.Leuven, Celestijnenlaan 200D, B-3001 Heverlee, Belgium 

Department of Physics and Astronomy, University of North Carolina, Chapel Hill, NC 27599-3255, USA 
Department of Astronomy & Astrophysics, Eberly College of Science, The Pennsylvania State University, 525 Davey 
Lab, University Park, PA 16802, USA 

Australian Astronomical Observatory, PO Box 296, Epping, NSW, 1710, Australia 

Department of Astrophysics, Faculty of Science, Radboud University Nijmegen, P.O. Box 9010, 6500 GL Nijmegen, 
NE 

Department of Physics, University of Warwick, Conventry CV4 7AL, UK 
9 Centre of Astrophysics Research, University of Hertfordshire, College Lane, Hatfield ALIO 9AB, UK 



S3 

(N 



4 



fjy . Argelander-Institut fur Astronomie, Auf dem Hiigel 71, 53121 Bonn, Germany 

O j The dates of receipt and acceptance should be inserted later 

6 

Key words binaries: spectroscopic - subdwarfs 

The project Massive Unseen Companions to Hot Faint Underluminous Stars from SDSS (MUCHFUSS) aims at finding 
hot subdwarf stars with massive compact companions (white dwarfs with masses M > 1.0 Mq, neutron stars or black 
holes). The existence of such systems is predicted by binary evolution calculations and some candidate systems have been 
found. We identified ~ 1100 hot subdwarf stars from the Sloan Digital Sky Survey (SDSS). Stars with high velocities 
have been reobserved and individual SDSS spectra have been analysed. About 70 radial velocity variable subdwarfs have 
■^j- | been selected as good candidates for follow-up time resolved spectroscopy to derive orbital parameters and photometric 

follow-up to search for features like eclipses in the light curves. Up to now we found nine close binary sdBs with short 
■ orbital periods ranging from ~ 0.07 d to 1.5 d. Two of them are eclipsing binaries with companions that are most likely 

\Q " of substellar nature. 
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1 Introduction Subdwarf binaries with massive WD companions may 

be candidates for supernova type la (SN la) progenitors be- 

| Hot subdwarf stars (sdO/B) are core-helium burning stars cause these systems lose angular momentum due to the emis-| 

^ located at the extreme blue end of the horizontal branch (see sion of gravitational waves and shrink. Mass transfer or the 

C$ Heber 2009 for a review). These stars with masses of about subsequent merger of the system may cause the WD to reach| 

half solar consist almost entirely of helium surrounded by a the Chandrasekhar limit and explode as a SN la. One of 

very thin hydrogen envelope only. The cause of the extreme the best known candidate systems for the double degenerate 

mass-loss in the red-giant phase necessary to form hot sub- merger scenario is the sdB+WD binary KPD 1930+2752 

dwarfs remains unclear. (Maxted et al. 120001 Geier et al. 120071 1 . 

A significant fraction (~ 50%) ofthe sdBs stars are in ^ ^ ^ d 201Qa | 2Q10bt ana lysed high resolution spec-| 

short period binaries (Maxted et al.[200l ; Napiwotzki et al. ^ Qf single . lined sdB binaries . Because the inc ii nations 

2004) with periods ranging from only 0.07 d to more than f , , • f 

1 ' & & j G f these systems are not known, additional information is 

10 d. These close binary sdBs are most likely formed by ■ ■ . ■ • • . , ^ ■ . 

3 . . 3 needed to derive companion masses. Accordingly, Geier et 

common envelope (CE) ejection (Han et al.| 2002|2003» . Be- d d 2010a [ | 2010bt per f orm ed a quantitative spectral analy- 

cause most of them are single-lined, only lower mass limits • , , . ■ , e ... , ■ . , . , 

& 3 sis and determined surface gravities and projected rotational 

have been derived from the binary mass functions consistent velocities Assuming synchronised orbits the masses and the 
with late main sequence stars of spectral type M or compact namre of ^ unseen companions were constrained . Surpris- 
objects like white dwarf s (WDs). ^ ^ masses of SQme companions are only consistent 

* Corresponding author: e-mail: geier@sternwarte.uni-erlangen.de with either massive white dwarfs (M > 1.0 M ), neutron 

stars (NS) or stellar mass black holes (BH). However, the 
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assumption of orbital synchronisation in close sdB binaries 
was shown to be not always justified and the analysis suf- 
fers from selection effects (Geier et al. 2010b). On the other 
hand, the existence of sdB+NS/BH systems is predicted by 
binary evolution theory (Podsiadlowski et al. 120021 Pfahl et 
al. 125531 Yungelson & Tutukov|2555j Nelemans|25l5J. The 
formation channel includes two phases of unstable mass 
transfer and one supernova explosion, while the fraction of 
systems formed in this way is predicted to be about 1 — 2% 
of all sdBs. 

If the companion were a neutron star, it could be de- 
tectable by radio observations as a pulsar. Coenen et al. 
(2011) searched for pulsed radio emission at the positions 
of four candidate systems from Geier et al. (2010b) using 
the Green Bank radio telescope, but did not detect any sig- 
nals. Most recently, Mereghetti et al. (1201 11 ) searched for X- 
ray signatures of mass transfer driven by weak stellar winds 
from the sdBs. Using the XRT instrument on board of the 
SWIFT satellite and targeting twelve binaries from the sam- 
ple of Geier et al. (2010b ), Mereghetti et al. (120111) did not 
detect any X-ray emission. 

In order to find sdBs with compact companions like mas-| 
sive white dwarfs (M > l.OM ), neutron stars or black 
holes we started a radial velocity (RV) survey (Massive Un- 
seen Companions to Hot Faint Underluminous Stars from 
SDSS0, MUCHFUSS, Geier et al. 1201 lall20TTbb . 

The same selection criteria that we applied to find such 
binaries are also well suited to single out hot subdwarf stars 
with constant high radial velocities in the Galactic halo like 
extreme population II and hypervelocity stars and led to a 
very interested spin-off project (Tillich et al. 1201 11 ). 



2 Target selection 



2.1 Colour and RV selection 



The target selection is optimised to find close massive com- 
pact companions to sdB stars. The SDSS spectroscopic database! The sample of promising targets consists of 69 objects 



2.2 Survey for RV variable stars 

Second epoch medium resolution spectroscopy (R = 1800— | 
4000) was obtained for 88 stars using ESO-VLT/FORS1, 
WHT/ISIS, CAHA-3.5m/TWIN and ESO-NTT/EFOSC2. Sec-| 
ond epoch observations by SDSS have been used as well. 
We discovered 46 RV variable systems in this way. 

The SDSS spectra are co-added from at least three indi- 
vidual "sub-spectra" with typical exposure times of 15 min, 
which are normally taken consecutively. Hence, those spec- 
tra can be used to probe for radial velocity variations on 
short timescales. We have obtained the sub-spectra for all 
sdBs brighter than g = 18.5 mag and discovered 81 new 
sdB binaries with radial velocity variations on short time 
scales (~ 0.03 d) in this way. In total we found 127 new RV 
variable hot subdwarf stars (see Fig. [I]). 

In addition, 20 helium-rich sdOs (He-sdOs) show RV 
variability. This fraction was unexpectly high since in the 
SPY sample only 4% of these stars turned out to be RV vari- 
able (Napiwotzki et al. 2008 ). However, it is not yet clear 
what causes this RV variability, since the orbital parameters 
of any such object couldn't be derived yet. This would be 
necessary to prove, that these He-sdOs are in close binaries. 

2.3 Selection of candidates with massive companions 

Numerical simulations were carried out to select the most 
promising targets for follow-up and estimate the probability 
for a subdwarf binary with known RV shift to host a massive 
compact companion. We created a mock sample of sdBs 
with a close binary fraction of 50 % and adopted the dis- 
tribution of orbital periods of the known sdB binaries. Two 
RVs were taken from the model RV curves at random times 
and the RV difference was calculated for each of the 10 6 bi- 
naries in the simulation sample. To account for the fact that 
the individual SDSS spectra were taken within short times- 
pans, another simulation was carried out, where the first RV 
was taken at a random time, but the second one just 0.03 d 
later (for details see Geier et al. 1201 lab . 



(Data Release 6) is the starting point for our survey. SdO/B 
candidates were selected by applying a colour cut to SDSS 
photometry. All point source spectra within the colours u — 
g < 0.4 and g — r < 0.1 were downloaded from the SDSS 
Data Archive Served About 10 000 hot stars were clas- 
sified by visual inspection. The sample contains 1100 hot 
subdwarfs (for details see Geier et al. 1201 lab . 

SdBs with radial velocities (RVs) lower than ±100 km s~ 
have been excluded to filter out such binaries with normal 
disc kinematics, by far the majority of the sample. Another 
selection criterion is the brightness of the stars since the 
quality of the spectra is not sufficient for faint stars. Be- 
cause of that most objects much fainter than g = 19 mag 
have been excluded. 



in total. These objects either show significant RV shifts (> 
SOkms- 1 ) within 0.03 d (52 stars) or high RV shifts (100- 
300 kms -1 ) within more than one day (17 stars). An exten- 
sion of our target selection to SDSS Data Release 7 is in 
progress. 



l|3 Sample statistics 



1 Sloan Digital Sky Survey 

2 das.sdss.org 



The classification of the hot subdwarf sample is based on 
existence, width, and depth of helium and hydrogen ab- 
sorption lines as well as the flux distribution between 4000 
and 6000 A. Subdwarf B stars show broadened hydrogen 
Balmer and He I lines, sdOB stars Hell lines in addition, 
while the spectra of sdO stars are dominated by weak Balmer| 
and strong He II lines depending on the He abundance. A 
flux excess in the red compared to the reference spectrum 
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Fig. 1 Highest radial velocity shift between individual spectra ( ARV) plotted against time difference between the corre- 
sponding observing epochs (AT). The dashed horizontal line marks the selection criterion ARV > 100 km s _1 , the dotted 
vertical line the selection criterion AT < 0.1 d. All objects fulfilling at least one of these criteria lie outside the shaded area 
and belong to the top candidate list for the follow-up campaign. The filled diamonds mark sdBs, while the open squares 
mark He-sdOs (Geier et al. (1201 lat . 



as well as the presence of spectral features such as the Mg I 
triplet at 5170 A or the Ca II triplet at 8650 A are taken as 
indications of a late type companion. 

In total we found 1100 hot subdwarfs. 725 belong to 
the class of single-lined sdBs and sdOBs. Features indica- 
tive of a cool companion were found for 89 of the sdBs and 
sdOBs. 9 sdOs have main sequence companions, while 198 
sdOs, most of which show helium enrichment, are single- 
lined (Geier et al. l2011al . 

4 Spectroscopy follow-up 

Follow-up medium resolution (R = 1200 — 4000) spec- 
tra were taken during dedicated follow-up runs with ESO- 
NTT/EFOSC2, WHT/ISIS, CAHA-3.5m/TWIN, INT/IDS, 
SOAR/Goodman and Gemini-N/GMOS. Orbital parameters 
of eight sdB binaries discovered in the course of the MUCH- 



FUSS project have been determined so far Geier et al. (1201 1 bl | 
l20TTcb . 

Because we deal with single-lined spectroscopic bina- 
ries, only their mass functions f m = M^ omp sin 3 i / (M comp +| 
M sdB ) 2 = PK z /2ttG can be calculated. Although the RV 
semi-amplitude K and the period P can be derived from 
the RV curve, the sdB mass M s dB> the companion mass 
A^comp and the inclination angle i remain free parameters. 
Adopting the canonical mass for core helium-burning stars 
MjdB = 0.47 Mq and i < 90° we derive a lower limit for 
the companion mass. 

With this minimum mass a qualitative classification of 
the companions' nature is possible in certain cases. For mini-| 
mum companion masses lower than 0.45 Mq a main se- 
quence companion can not be excluded because its lumi- 
nosity would be too low to be detectable in the optical spec- 
tra (Lisker et al. |20051 >. The companion could therefore be 
a compact object like a WD or a late main sequence star. 
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If the minimum companion mass exceeds 0.45 Mq and no 
spectral signatures of the companion are visible, it has to 
be a compact object. If the mass limit exceeds 1.00 Mq or 
even the Chandrasekhar limit (1.40 Mq) the existence of a 
massive WD or even an NS or BH companion is proven. 

The derived minimum companion masses of seven bi- 
naries from our sample are similar (0.32 — 0.41 Mq). From 
these minimum masses alone the nature of the companions 
cannot be constrained unambiguously. However, the fact thafl 
all seven objects belong to the sdB binary population with 
the highest minimum masses illustrates that our target se- 
lection is efficient and singles out sdB binaries with massive 
companions (see Geier et al. 1201 lbb . 

5 Photometry follow-up 

Photometric follow-up helps to clarify the nature of the com-| 
panions. Short period sdB binaries with late main sequence 
or substellar companions show variations in their light curves| 
caused by the irradiated surfaces of the cool companions 
facing the hot subdwarf stars. If this so-called reflection 
effect is present, the companion is most likely a main se- 
quence star. If not, the companion is most likely a compact 
object. In the case of the short period system Jl 138— 0035 a 
light curve taken by the SuperWASP project (Pollacco et al. 
2006) shows no variation exceeding ~ 1%. The companion 
is therefore most likely a white dwarf (Geier et al. 1201 lbl l. 

We obtained follow-up photometry with the Mercator 
telescope and the BUSCA instrument mounted on the CAHA-| 
2.2m telescope. In this way we discovered the first eclipsing 
sdB binary J0820+0008 to host a brown dwarf companion 
with a mass ranging from 0.045 to 0.068 Mq (Fig. [2] Geier 
et al. !2011cl >. 

A very similar eclipsing system (J 1622+4730) was dis- 
covered serendipituously (see Fig. [3]). A preliminary analy- 
sis shows that the orbital period is very short (~ 0.07 d) and 
the RV semi-amplitude quite low (~ 47kms~ 1 ). The com- 
panion is most likely a substellar object as well. The high 
success rate in finding these objects shows that our target 
selection not only singles out sdB binaries with high RV- 
amplitudes, but also systems with very short orbital periods 
and moderate RV-amplitudes. 

Most recently, we detected p-mode pulsations in the sdB 
JO 120+3950 (FBS 01 17+396, Geier et al. 120 11 al l as well as 
a longer trend indicative of a reflection effect in a light curve 
taken with BUSCA. Only a few of the known short-period 
sdB pulsators (sdBV r ) are in close binary systems. 

6 Substellar companions 

The existence of eclipsing sdB+ late dM binaries of HW Vir 
type with short orbital periods (0.1 — 0.26 d) and compan- 
ion masses between 0.1 Mq and 0.2 Mq (For et al. 2010; 
0stensen et al. 2010 and references therein) shows that stars 
close to the nuclear burning limit of ~ 0.08 Mq are able to 



help eject a common envelope and form a hot subdwarf. 
Substellar companions to sdB stars have been found using 
the light travel time technique (Schuh 2010 and references 
therein). However, these systems have wide orbits and none 
of these companions influenced the evolution of its host star. 

In the course of the MUCHFUSS project we discovered 
two sdBs with most likely substellar companions in close 
orbits. These companions evidently interacted with the sdB 
progenitor stars and caused the ejection of the common en- 
velopes. Some theoretical models indeed predict such an in- 
teraction between planets or brown dwarf companions and 
their nearby host stars to be a possible formation channel 
for hot subdwarfs and helium white dwarfs (Soker 1998; 
Nelemans & Tauris 1998). Our finding can be used to con- 
strain such models and learn more about the role substellar 
companions play in the formation of single and close binary 
sdBs. 

7 Summary 

The MUCHFUSS project aims at finding hot subdwarf stars 
with massive compact companions. We identified 1100 hot 
subdwarfs by colour selection and visual inspection of the 
SDSS-DR6 spectra. The best candidates for massive com- 
pact companions are followed up with time resolved medium] 
resolution spectroscopy. Up to now orbital solutions have 
been found for eight single-lined binaries. Seven of them 
have large minimum companion masses compared to the 
sample of known close binaries, which shows that our target 
selection works quite well. However, it turns out that our se- 
lection strategy also allows us to detect low-mass compan- 
ions to sdBs in very close orbits. We discovered an eclips- 
ing sdB with a brown dwarf companion and a very similar 
candidate system in the course of our photometric follow- 
up campaign. These early results encourage us to go on, 
because they demonstrate that MUCHFUSS will find both 
the most massive and the least massive companions to sdB 
stars. 
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Fig. 2 Phased R-band light curve of J0820+0008 taken with the Mercator telescope. A light curve model is fitted to the 
data and residuals are given below. Primary and secondary eclipses can be clearly seen as well as the sinusoidal shape 
caused by the reflection effect (Geier et al. l201 let . 
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Fig. 3 Phased light curves of J1622+4730 taken with BUSCA (UV,B,R,IR-band). Although the eclipses are not total, the 
light curve looks very similar to the one of J0820+0008 (see Fig.0. 
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